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Abstract 

We discuss exclusive elastic double diffractive axial- vector Xc(l + ) meson production in proton- 
antiproton collisions at the Tevatron. The amplitude for the process is derived within the kf 
factorisation approach with unintegrated gluon distribution functions (UGDFs). We show that 
the famous Landau- Yang theorem is not applicable in the case of off-shell gluons. Differential 
cross sections for different UGDFs are calculated. We compare exclusive production of Xc(l + ) and 
Xc(0 + )- The contribution of x c (l + ) to the J/^+j channel is smaller than that of the Xc(0 + ) decay, 
but not negligible and can be measured. The numerical value of the ratio of the both contributions 
is almost independent of UGDFs modeling. 
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I. INTRODUCTION 



The central exclusive production of mesons has been recently revived. This is essentially 
because of two reasons. 

Firstly, the theoretical QCD inspired approach has been developed. This is because of 
interest in the double diffractive production of the Higgs boson firstly proposed by A. B. 



, 0, i] (KKMR) as an alternative 



Kaidalov, V. A. Khoze, A. D. Martin and M. G. Ryskin 
to inclusiveproduction for Higgs searches. In principle, very similar methods can be used 
for scalar [7|, pseudoscalar 



axial-vector and tensor mesons. The situation for vector 
meson production is somewhat different. Here the dominant mechanism is photon-pomeron 
(pomeron- photon) fusion j^] or pomeron-odderon (odderon-pomeron) fusion 03- Recently a 
/^-factorization approach has been used to calculate exclusive T production [11] at Tevatron. 
Secondly, some experimental efforts have been done to facilitate real measurements at 



Tevatron 12| and in the future at RHIC_ 13|_and LHC 14 
Tevatron have been presented recently 



12, 15 



Some preliminary results from 




FIG. 1: The sketch of the bare QCD mechanism for diffractive production of the Xc(l + ) meson. 
Some kinematical variables are shown in addition. 

In the present paper we are concentrated on the exclusive double-diffractive production 
of axial- vector x c (l + ) mesons 1 . Here the dominant mechanism is the two-pomeron fusion, 
which in the QCD language is a "fusion" of two QCD ladders. The mechanism is shown 
schematically in Fig. [TJ Compared to the Higgs production, where a hard scale is guaranteed 
by the large mass of the Higgs boson, here the natural scale (mass of the Xc(l + ) meson) 
is much lower and the method proposed by KKMR is a bit questionable as a big part 
of the strength may come from the region of relatively small gluon transverse momenta. A 
pragmatic solution is to use nonperturbative models of UGDFs instead of the pQCD inspired 
KKMR procedure (for more details see 0]). 

The situation with the axial-vector production is new compared to both zero-spin case 
(scalar 0], pseudoscalar [s| mesons) as well as to the vector meson production where the 



Some general aspects of this process were discussed previously in Refs. 
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vector meson is dominantly transversely polarized 0, [ll|, at least, for small transferred 
four-momenta in the nucleon lines. The axial-vector meson, as it will be discussed here, can 
be polarized both transversely and longitudinally. We shall calculate the cross section for 
different polarisation states of the Xc(l + ) meson. 

There is interesting theoretical aspect of the double diffractive production of the Xc(l + ) 
meson. The coupling g*g*xd^ + ) (see Fig. [T]) vanishes for on-shell gluons (so-called Landau- 
Yang theorem). According to the original Landau- Yang theorem [16j the symmetries under 
space rotation and inversion forbid the decay of the spin-1 particle into two (on-shell) spin-1 
particles (two photons, two gluons). The same is true for the fusion of two on-shell gluons. 
The symmetry arguments cannot be strictly applied for off-shell gluo ns. This fact has been 



already explored in inclusive production of Xc(l + ) [13> UM, Il9|, l2Q|, in the production of 
spin-1 glueballs j2l|, and recently in the studies of decays of hypothetical Z' bosons into 
pair of standard Z bosons [221 ] . One of the goals of our paper is to confirm explicitly that 
the Landau- Yang theorem is violated by virtual effects in diffractive production of Xc(l + ) 
leading to very important observational consequences. In our approach the off-shell effects 
are treated explicitly. For comparison, in the standard KKMR approach the corresponding 
cross section would vanish due to their on-shell approximation. The measurement of the 
cross section can be therefore a good test of the off-shell effects and, consequently, UGDFs 
used in the calculation. 

At the Tevatron the Xc mesons are measured through the 7 + J/\I/ decay channel. The 
axial-vector % c (l + ) meson has a large branching fraction of the radiative decay % c (l + ) — > 



7 + J/ip (BR = 0.36 [23j])- This is much bigger than for the scalar Xc(0 + ) where it is 



only about 1 % [23J. Therefore, the discussed off-shell effects may be very important to 
understand the situation in the 7 + J channel observed experimentally. 



II. FORMALISM 

The kinematics of the process was already discussed in our previous paper on Xc(0 + ) Q. 
Here we discuss only details of the matrix element for exclusive Xc(l + ) production. This is 
derived for the first time, at least in the generalized KKMR approach. 



A. General \cj production amplitude 

In the following we employ the general Kaidalov-Khoze-Martin-Ryskin approach [H, 0, 0] , 
and write the amplitude of the exclusive double diffractive color singlet production pp 

ppXcj as 

= S - • TT^j^a / d\ , t V^( qu q 2 ,p M ) 
fgfi( x i> x ii <2o,t> Qi,v ^i)/s,2 / ( x 2, x' 2 , ql t , q\ p t 2 ) 

X 2 2 2 • v > 

%,t Ti, t 

In this expression fg^/ixi^x'^ql^qf^ti) are off-diagonal unintegrated gluon distributions. 
In the general case we do not know the off-diagonal UGDFs very well. In Ref. 0, @] we 
have proposed a prescription how to calculate the off-diagonal UGDFs with the help of their 
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diagonal counterparts: 



f?J = \Jff\^qit^l)-ti 2 \^qt t ^ 2 ) ■ F 1 (t 2 ) , (2.2) 

where and Fife) are the isoscalar nucleon form factors. In the present work we 

shall use a few sets of unintegrated gluon distribution functions (UGDFs), which aim at 
description of phenomena where small gluon transverse momenta are involved. Some details 
concerning these distributions can be found in Ref. 24|. We shall follow notations there. 

Following our previous work Q, the vertex factor Vj lC2 = Vj lC2 (qf t , g| t , Pl I)t ) in Eq. (12. ip 
describing the coupling of two virtual gluons to Xcj- m eson follows from 

Vr 2 (QhQ2) =V{qq^ X cj)**T(qi,q2), (2.3) 

where V(qq — > Xcj) is the operator that projects the qq pair onto the charmonium bound 
state (see below), \I/ ClC2 (gi, q 2 ) is the production amplitude of a pair of massive quark q and 
antiquark q with momenta k±, k 2 , respectively. 

Within the quasi-multi-Regge-kinematics (QMRK) approach [25| we have 

*( Cl , c 2 ; i, k; qi , q 2 ) = -g\tt]t%b{k x , k 2 ) - t c k 2 tfb(k 2 , k x )), a s = (2.4) 

where t c are the colour group generators in the fundamental representation, b, b are the 
effective vertices arising from the Feynman rules in QMRK 

b{ki, k 2 ) = 7 7 TTyi 2'' °(«i> = 7 7 rr^ 2^ ■ 2 - 5 

While projecting on the color singlet the ggg-vertex contributions disappear from the result- 
ing matrix element, so we did not write them explicitly in Eq. (12. 5p . Taking into account 
standard definitions of the light-cone vectors n + = p 2 /E cms , n~ = pi/E cms and momentum 
decompositions q\ = x\pi + qx st , q 2 = x 2 p 2 + q 2jt and using the gauge invariance property 
(Gribov's trick) one gets the following projection 

^T/Cic 2 _ „MT/cic 2 _ n 
9l V J,au — V2 V J,au — U > 



J, \av - Hi > J, fiv 

1 ,,/' ; 



4 tf t 



Vr 2 ( qi ,q 2 ) = n;n-V^( qi ,q 2 ) = -^^V^( qi , q 2 ). (2.6) 

Since we adopt here the definition of the polarization vectors proportional to gluon transverse 
momenta 91/2,$, then we must take into account the longitudinal momenta in the numerators 
of vertices (12.51) . 

Projection of the hard amplitude onto the singlet charmonium bound state V^ C2 is given 
by the 4-dimensional integral over relative momentum of quark and antiquark q = (ki — k 2 )/2 

.0 q 2 



V^(qu 02) = V(qq -> X cj) • = ^ ' E E ^ / (^fl* M 

i,k L Z ,S Z 

x $ l =i,lM '(L = l, L z ; 5=1, S Z \J, J z ) (3z, 3/c|l) Tr {* t %V>s=i A } , (2.7) 

9 2^[Vi* |> (gi _ A , i)2 _ m2 7Mj > Vii \^( gi _ A;2 )2_ m 2^ 
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where the function $L=i i L z (q) is the momentum space wave function of charmonium, and 
for a small relative momentum q the projection operator Vs=i,s z has the form 



Vs=i A = ^{k - m)^-(k + m) . (2.8) 

Since P-wave function *&l=i,l„ vanishes at the origin, we may expand the trace in Eq. (I2.7P 
in Taylor series around q = 0, and only the linear terms in q a in the trace survive. This 
yields an expression proportional to 

/ Wf qa ^ L ^ = -^^^ ^ 

with the derivative of the P-wave radial wave function at the origin 72/(0) whose numerical 
values can be found in Ref. 27] . The general P-wave result (12.71) may be further reduced 
by employing the Clebsch-Gordan identity which for the vector Xcj=i charmonium reads 

r/4 = ^(l.L.ji.ftii, j t )^(L M y(s M )=-i^l^^ep(j t ) . (2.io) 



L.,S. 



B. gg — > Xc(l + )-vertex function 

Summarizing all ingredients above in Eqns. (12.61) . (12.71) . (12. 9p and (12.101) . we get the vertex 
factor in the following covariant form 
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qlAMPi - ^M, t + qlt) - (2.ii) 
P1P2 ( ^(2?2, t (gi? 2 ) - (qi,tg2,t)(qlt + qlt)) - qtt^qlM^) - (qi,tq2,t)(qlt + qlt)) 



The general vertex function (12.111) possesses the Bose symmetry under simultaneous per- 
mutation of gluon momenta q\ <-> q2 and polarisation vectors n + <-> n~ defined in Eq. (12.61) . 
or, equivalently, under simultaneous permutations of protons (pi <-> P2) and gluons (both 
transverse q\j <->• q2,t and longitudinal X\p\ <-» X2P2) momenta. 

We write the decomposition of the polarisation vector of a heavy meson with a given 
helicity A = 0, ±1 as 

e P{P,\) = {l-\\\)4-±={\4 + i\\\4), < = <n^=g a p, e M (A)e*(A') = —5 XX 

In the c.m.s. frame we choose the basis with collinear 113 and P vectors (so, we have 
P = (P, 0, 0, P z ), P z = |P| > 0) as a simplest one 

n? = (0,1,0,0), 4 = (0, 0, 1, 0), 4 = i(|P|,0,0,P), |P| = y/E*-M\ (2.12) 
Note, that we choose n 2 to be transverse to the c.m.s beam axis (see Fig. [2]), while n 1; n 3 
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FIG. 2: Coordinate basis in the center-of-mass system of incoming protons pi,2- 



are turned around by the polar angle ip = [0 ... 7r] between P and the c.m.s. beam axis. In 
the considered basis {rii, n 2 , n 3 } we have the following coordinates of the incoming protons 

fs fs 
Pi = — (1, -sin-0, 0, cos-0), p 2 = —(1, siiaip, 0, -costp) . (2.13) 

The gluon transverse momenta with respect to the c.m.s. beam axis are 

qi,t = (0, Qf it cosip, Qt, Qltsimfj), q %t = (0, Q^cos^, -Q\, Q^ t sin^), 

where Q\i 2v ±Qt are the components of the gluon transverse momenta in the basis with 
the z-axis collinear to the c.m.s. beam axis. 

From definition (12.131) it follows that energy of the meson and polar angle ip are related 
to covariant scalar products in the considered coordinate system as 

B= (ga+fafl, cosi , = (P' p '>'(rf\ ^ =(»*>.) -(*"■), ( 2.i4) 

y/S V S I P I V S 

Further, we also see that from q 1 = X\P\ + qi tt , q 2 = x 2 p 2 + q 2 ,t and qi + q 2 = P we have 

E + |P| cosip E — |P| cosip 
xi = -= , x 2 = y= . (2.15) 



Relations (12.141) and (12.151) show that the interchange of proton momenta p\ ^ p 2 is equiv- 
alent to the interchange of the angle tp «-> tp ± 7r, i.e. simp «-> — sin-^ and cosip <-> — cos^, 
simultaneously. The last permutation also provides the interchange of the longitudinal com- 
ponents of gluons momenta X\ «-> x 2 . 

Conservation laws provide us with the following relations between components of gluon 
transverse momenta and covariant scalar products 



qi t + (?i,t52,t) ^ ql t + («i,t?2,t) ^ \Jqltqlt-(qi,tq2,? 2 



Qi,t — ~ ip| . ^ ; = .p. / ? Qr = ^ jp-j si g n(gr 

|P| sin?/; |P| sm^ | Pi | 

P 2 = -|P t | 2 = -|P| 2 sin 2 V = ql t + ql t + 2{q 1>t q 2>t ), qj /2it = -|qi/ 2 , t | 2 , 



where |P 4 | = |P||sin^| is the meson transverse momentum with respect to z-axis. The 
appearance of the factor sign(Q^) guarantees the applicability of Eq. (12 . 1 6[) for positive and 



negative Q\. Note that under permutations qi tt <-> g 2 ,t implied by the Bose statistics the 
components interchange as Qf t <-> Qf t and «-> — <5t- In our notations the quantity sinip 
plays a role of the noncollinearity of meson in considered coordinates. A straightforward 
calculation leads to the following vertex function in these coordinates 



V J=l,\ 



2 rciC2 



G 



ft'(O) 



MnN c \P t \(M* - ql, - ql,Y\^2 



- 92, t )(M2,t)sign(sinV0 



+ Kilt + ?2,t) I [qi,* x <te,t] x ni| sign(Qf) sign(cos^) 



+ 



(2.16) 



+ (1 - W)(sL + ?|,t) |[qi,t x <j2,t] x n a| sign(g^) sign(sin^) 



where 



|[qi,t x q 2 , t ] x nil = \ q\ t q\ t - {q\,tq%t) 2 | cos^|, 



|[qi,t x q 2 ,t] x n 3 | = —\/ql t ql t - (qi,tq2,t) 2 | sin^|. 



The amplitude (12.161) explicitly obeys the Bose symmetry under the interchange of gluon 
momenta and polarisations due to resulting simultaneous permutations cosip «-> —cosip, 
sinip «-> — sin^ and ~Qt 2 - 

A short inspection of Eq. (12.161) shows that 



(2.17) 



when gi t — > or q^ t — > 0. It shows that gluon transverse momenta (gluon virtualities) are 
necessary to get a nonzero cross section. It also means that the amplitude and the cross 
section are sensitive to larger values of gluon transverse momenta than e.g. in the case of 
X c (0 + ) production. 

It follows from the conservation laws that 

Qu + Pit = -lot, q» + Pat = lot, Pt = -(p'u + Pat) 
Let us consider firstly the limit of the "coherent" scattering protons p' lt = p' 2t = Pt, so 

Qit = -(p* + <?ot), qn = -{pt - 9ot), P t = —2pt, P v t =0- (2.18) 
The production vertex (I2.16P in this limit and considered coordinates has a form 

mIn c w - Wh ti t )A mqlt - p2 ^ si ^ sin ^ + 



V^lMtAt,Pt) = -i6g 2 S ClC2 

+ (Pt+iMt 



A cos ip H — ]^ (1 — | A|) sin -0 



(2.19) 



2 We are thankful to M. G. Ryskin for enlightening correspondence on the issues of Bose symmetry of 
production amplitude. 
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This vertex is antisymmetric w.r.t. simultaneous interchanges q$ t — q$ t and q$ t <-> — g w : 

VjllMt, 5m, Pt) = "^S!a(-9&. (2-20) 
In the considered "coherent" limit (12.181) the integrand of the diffractive amplitude 

Vj=lx(<fon qot,Pt) ■ fg!i(xi,x'i, qju (pt + got) 2 , t 1 )f° f 2 f {x 2 , 4, g g if , (p t - g 0f ) 2 , t 2 ) 

will be antisymmetric only if x\ = x 2 = E/y/s = x (while x[ ~ x 2 <C £1,2)) i>e. in the 
case when y = 0, while the deviation from zero at y 7^ manifests the violation of Regge 
factorization which was used to examine this limit earlier [jjj]. So, the diffractive amplitude 
in this case 

2 „2 \ f(„ Jl Jl 



KA T? (+ \T? U \ f A x A V V J=l(<ftt> QouPt) " f(?, 9o,t> ?l,t)/ ( X > %,t> %t) n 

-M^o ~ F^F^h) / dg ^4 -2-2-^ ! ! = 

J ( iOt < ilt t i2t 



In the forward limit p t — ► (which is the particular case of coherent one) the amplitude 
turns to zero at any y. Indeed, we have P t — > and sin ip —>■ ±0 and the amplitude turns 
into 



VjL c UQot,Q y ovPt - 0) = -16^7^(0),/— . x (2 . 2 1) 



x ( M 2 !°2 g 2 si | j l A |goi si g n ( sin ^)U-^o,7r + Aq&sign(cosV>)U-o,7r J 



As in the previous case, it is obviously antisymmetric under interchanges q$ t <-> —q^ t and 
q^ t <-> — got- Since in this case qu = —qot, qit = qot, then the diffractive amplitude has an 
antisymmetric integrand and turns to zero 

M pt ^ ~ Fi(ti)Fi(t 2 ) / dcfadqlt g ; ; ; = 0. 

J %t 

This explicitly confirms the observation made in Refs. 0, [H| 3 . 

It is also possible to express the results in terms of transverse 3-momenta of fusing off- 
shell gluons I qi,* I and |q2,t|, and the angle between them in the center-of-mass system of 
colliding nucleons with z-axis fixed along meson momentum P. In this case, summing the 
squared matrix elements over meson polarizations we get the expression 



In I 2 

iqi,*i 


q2,*l 2 


^(|qMl 2 + 


qi,t 


2\ 2 • 

) sm 


2 + M 2 (|q M 


2 + 


02,* 


2 "2|q M | 


|q 2)t | cos0) 


(|q M 


2 + 


q 2 , t | 2 + M 2 ) 4 



equal (up to different normalisations of gluon polarization vectors) to the one derived in 
Ref. El. 



We are grateful to V. A. Khoze for very interesting and helpful correspondence on this problem. 
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C. Three-body phase space 

At high energies and small momentum transfers the phase space volume element can be 



written as 28 



d3pS = ^dhdhd^d&d® 5 (s(l - £0(1 - 6) " M 2 ) , (2.22) 

where £i, £2 are longitudinal momentum fractions carried by outgoing protons with respect to 
their parent protons and the relative angle between outgoing protons $ 6 (0, 2ir). Changing 
variables (£1,^2) — > (xf,M 2 ) one gets 

d 3 PS = -J—dtjdh—, dXF rf$ . (2.23) 



III. RESULTS 

Let us start from presenting the differential cross sections. In Fig.[3]we show distributions 
in rapidity y for different UGDFs from the literature. The results for different UGDFs signif- 
icantly vary. The biggest cross section is obtained with BFKL UGDF and the smallest one 
with Gaussian UGDFs. The big spread of the results is due to quite different distributions 
of UGDFs in gluon transverse momenta qu, q2t, although when integrated over transverse 
momenta distributions in longitudinal momentum fractions Xi,X2 are fairly similar. 




y y 



FIG. 3: Distributions in rapidity of Xc(l + ) meson (left panel) and Xc(0 + ) meson (right panel) 
for different UGDFs. Dash-dotted line corresponds to BFKL UGDF, long-dashed line - GBW, 
short-dashed line - KL, and two solid lines - Gaussian UGDFs for o$ = 0.5 GeV 2 (upper line) and 
(To = 1.0 GeV 2 (lower line). 

Comparing the left and right panels, the cross section for the axial-vector x c (l + ) produc- 
tion is much smaller (more than an order of magnitude) than the cross section for the scalar 
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Xc(0 + ) production. This is related to the Landau- Yang theorem, which "causes" vanishing 
of the cross section for on-shell gluons. For axial-vector quarkonia the effect is purely of 
off-shell nature and is due to the interplay of the off-shell matrix element and off-diagonal 
UGDFs. This interplay causes a huge sensitivity of differential distributions to UGDFs 
observed in Fig. [3J 




FIG. 4: Distribution in tip of Xc(l + ) meson (left panel) and Xc(0) meson (right panel) for different 
UGDFs. 

In Fig. H] we show corresponding distributions in t — t\ or t — £2 (identical) again for 
different UGDFs. Except of normalisation the shapes are rather similar. This is because of 
the t\ and £2 dependencies of form factors, describing the off-diagonal effect, taken the same 
for different UGDFs. 

In Fig. [5] we show the correlation function in relative azimuthal angle between outgoing 
protons. The shapes of the distributions are almost independent of UGDFs. In the case when 
energy resolution is not enough to separate contributions form different states of Xc (Xc(0 + ), 
Xc(l + )> Xc(2 + )) the distribution in relative azimuthal angle may, at least in principle, be 
helpful. 

Summarizing differential distributions, the cross sections (especially their absolute nor- 
malisation) strongly depend on the model of UGDF. In spite of the huge uncertainty in 
predicting the absolute cross section it becomes obvious that the cross section for x c (0 + ) 
is much bigger than the cross section for the Xc(l + ) production. This result could be ex- 
pected based on the Landau- Yang theorem. However, the size of the suppression cannot be 
predicted without actual calculations. 

In Table 1 we have collected cross sections integrated in y, ti, t 2 , 4> over the full phase space 
for the Tevatron energy W = 1960 GeV. More than an order of magnitude suppression of 
X c (l + ) relative to Xc(0 + ) can be seen by comparing numbers in appropriate columns 4 . 



4 This ratio should be only weakly modified by absorption effects. 
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FIG. 5: Distribution in relative azimuthal angle <3? of Xc(l + ) (left panel) and Xc(0 + ) (right panel) 
meson production for different UGDFs. 



TABLE I: Integrated cross section a tot (in nb) for exclusive Xc(0 + ) and Xc(l + ) production for 
different UGDFs and the Tevatron energy W = 1960 GeV. Branching ratios of radiative decays 
were taken from [29]: BR(x c (0+) J/¥ 7) = 0.0128 and BR( Xc (l + ) -> J/* 7) = 0.36. 



UGDF 


Xc(0+) 


Xc(l + ) 


ratio 

BR-< Hot ( X cU + )) 
BR-a tot ( Xc (0+)) 


&tot 






BR-d tot 




KL 


55.2 


0.7 


0.5 


0.2 


0.3 




GBW 


160 


2 


4.2 


1.5 


0.8 




BFKL 


1200 


15.4 


14.2 


5.1 


0.3 




Gauss, 














= 0.5 GeV 


26 


0.3 


0.2 


0.09 


0.3 




Gauss, 












O"0 


= 1.0 GeV 


2.2 


0.03 


0.02 


0.006 


0.2 



The best method to measure Xc mesons at the Tevatron is via 7 + J/\I/ decay channel. 
All P-wave x c -quarkonia decay into this channel. However, the branching fractions to this 
channel are very different [23|. While the branching fraction for Xc(0 + ) is very small (of 
the order of 1 %), the branching fraction for % c (l + ) is one and half order of magnitude 
larger. In the third and fifth columns we present the total cross sections multiplied by 
the appropriate branching fractions. After multiplying the cross section by the branching 
fraction for 7 + J/ty decay the situation somewhat changes, i.e. now x c (l + ) becomes closer 
to Xc(0 + ) - BR-<7to((xc(l + )) is about two times smaller than that for Xc(0 + ) for all UGDFs 
used in our calculation, so it seems to be almost model independent statement. 
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In a preliminary analysis the CDF collaboration 15| assumes that the observed strength 
comes dominantly from Xc(0 + ), thus conforming results of our investigation. In order to 
make comparison with the experimental results one would still need to include experimental 
cuts on lepton and photon rapidities and transverse momenta. Also including absorption 
effects may be important as slightly larger absorption can be expected for x c (l + ) (harder 
distributions in t\ and t 2 - see Fig. @J. These points need further studies. 

Our calculation suggests that the inclusion of Xe(l + ) m the experimental analysis is 
not negligible and may be necessary. The present energy resolution does not allow for 
separating different P-wave quarkonia. Perhaps, looking to other decay channels may help 
in disentangling the contributions from different states and allowing for extracting the cross 
sections separately for each of them. 

Another interesting option to shed more light to the problem is to study the angular 
distributions of outgoing J /if) in the Xc rest frame. Different states should have, in principle, 
different distributions. We leave the analysis of those distributions for a separate study. 



TABLE II: Integrated cross section otot (in nb) for exclusive x c (l + ) production at different energies. 



UGDF 


RHIC 


Tevatron 


LHC 




KL 


0.05 


0.5 


1.7 




GBW 


0.04 


4.2 


73.1 




BFKL 


0.07 


14.2 


1064 




Gauss, 










= 0.5 GeV 


0.007 


0.2 


2.5 




Gauss, 








O"0 


= 1.0 GeV 


0.0005 


0.02 


0.2 



Finally in Table 2 we present the total cross sections for x c (l + ) also for RHIC and LHC 
energies. The question of separation of different Xc states should be similar, except that 
other decay channels should be available [Lj, uA ■ 



IV. CONCLUSIONS AND DISCUSSION 



Our results can be summarized as follows: 

We have derived the QCD amplitude for exclusive elastic double diffractive production 
of axial- vector x c (l + ) meson. According to the Landau- Yang theorem the amplitude van- 
ishes for the fusion of on-shell gluons. In the present analysis we have generalized the 
formalism proposed recently for diffractive production of the Higgs boson. We have derived 
corresponding g*g* — > x c (l + ) vertex function. Our effect is purely off-shell type, i.e. re- 
quires off-shell gluons, which demands nonvanishing transverse momenta of gluons in the 
high-energy regime. 

We have calculated the corresponding differential cross sections. Different unintegrated 
gluon distributions from the literature have been used. The absolute cross section is very 
sensitive to the choice of UGDF in contrast to the shapes of distributions. The predicted 
total (integrated over phase space) cross section, obtained from the bare amplitude, is from 
a fraction to several nanobarns, depending on the model of UGDFs. This is one and a 
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half order of magnitude less than a similar cross section for Xc(0 + ) [7]. This is a di- 
rect consequence of the Landau- Yang theorem. However, because the branching fraction 
BR(xc(l + ) — > J /i> + 7) 3> BR(xc(0 + ) — > J /if) + 7), one may expect a different situation in 
the J /if} + 7 channel. This has an analogy with the inclusive production of P-wave quarko- 
nia, where the signal (in the J /if> + 7 channel) of x c (l + ) is larger than that for x c (0 + ). We 
have observed that BR-er tot (x c (l + )) is (only) several times smaller than that for Xc(0 + ) for 
all UGDFs used in our calculation. 

Moreover, we have also calculated [30J differential cross sections for different spin polariza- 
tions of x c (l + ). The integrated cross section for spin polarization A = ±1 is approximately 
an order of magnitude greater than that for the A = polarization. Similar observation has 
already been made in Refs. 0,11] and verified by the WA102 data for /i(1285), /i(1420) pro- 
duction [3JJ. The ratio of the cross sections integrated over the phase space is only weekly 
dependent on UGDFs but strongly depends on t\ and t 2 . 

In the present analysis we have neglected the absorption effects. The latter clearly go 
beyond the scope of the present analysis. At the Tevatron energies they lead, however, to a 
large damping of the cross section. In zeroth approximation they can be taken into account 
by multiplying the cross section by a so-called soft survival probability [l], 0] . At the Tevatron 
energies the soft survival probability is of the order of 0.1 [3, 0]. A better approximation is 
to convolute the bare amplitude with nucleon-nucleon elastic (re)scattering amplitude (see 
e.g. Refs. (U [ll|, 32]). We leave the inclusion and discussion of the absorption/rescattering 
effects for a separate analysis, including all quarkonium states (Xc(0 + ), Xc(l + ), Xc(2 + )). 
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